Spin Mixing in Spinor Fermi Gases 
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We study a spinor fermionic system under the efTect of spin-exchange interaction. We focus on 
the interplay between the spin-exchange interaction and the effective quadratic Zeeman shift. We 
examine the static and the dynamic properties of both two- and many-body system. We find that 
the spin-exchange interaction induces coherent Rabi oscillation in the two-body system, but the 
oscillation is quickly damped when the system is extended to the many-body case. 
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I. INTRODUCTION 

Spinor quantum gas has received tremendous attention 
ever since the first creation of a spin-1 Bose-Einstein con- 
densate in an optical trap [l| ■ Most of these studies deal 
with spinor Bose gases. The most salient feature of spinor 
condensates is the presence of spin-exchange interaction 
which drives coheret spin-mixing dynamics. A natural 
question can be asked is: can similar behavior be ob- 
served in a quantum degenerate fermionic system? The 
current work is an attempt to address this question. 

Quantum degenerate Fermi gases have been realized in 
many laboratories 0411 ■ Different with the normal elec- 
tronic system, many fermionic atoms have spins higher 
than 1/2 in their lowest hyperflne manifold. These large- 
spin ultracold Fermi gases provide us a unique opportu- 
nity to investigate exotic many-body physics [6|-la[, and 
have stimulated a great deal of theoretical interest p3l-[l2| . 
Considerable experimental progress has also been made 
recently in the system of *^Sr (/ = 9/2, where / is total 
hyperfine spin) jll and ^^^yb (/ = 5/2) [3. Both ^^Sr 
and ^^'^Yb have an alkaline-earth-like atomic structure 
with all electron shells filled, thus their hyperfine spins 
completely come from rmclear spins. This will lead to 
a spin-independent atom-atom interaction since the nu- 
clear spins are deep within the atom. In these alkaline- 
earth-like fermionic system, the spin-independent inter- 
action can give rise to the so-called SU{N) symmetry, 

with N ^ 2f + 1 However, for a more general 

cold Fermi system including non-alkaline-earth atoms 
with large spins, the SU (N) symmetry may not be re- 
served. As a result, one can find a more rich phases 
diagram in its ground state [18]. In a seminal work pub- 
lished very recently [19 ], a spinor Fermi gas of ^"K was 
realized. By taking advantage of the spin conservation, 
the effective spin of the system can be tuned from 1/2 to 
9/2. 

In this article, we focus on the simplest large-hyperfine- 
spin systems with f — 3/2 with four internal components. 
This can be either true hyperfine spin (e.g., ^^^Cs, ^Be 
and ^°^Hg) or effetive spin such as realized in Ref. [l9| . 
If we consider atoms with nonzero electron spins due to 
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FIG. 1: (Color Online) Schematic of a spin-mixing process. 



partially filled electron shells, then the interaction among 
atoms will be spin-dependent. One of the key features in 
this kind of systems is that there will be spin-exchange 
interactions which constantly mix different spin compo- 
nents. For example, two atoms with respective hyperfine 
spins —1/2 and +1/2 interact and become two atoms 
with hyperfine spins —3/2 and -1-3/2, as schemcatically 
depicted in Fig.[T] A similar case of spin mixing has been 
well studied for spinor condensates |22H25j . 

The paper is organized as follows: We present the 
model Hamiltonian in Sec. II. The ground state proper- 
ties and the spin- mixing dynamics are discussed in Sees. 
Ill and IV, respectively. Finally, we conclude in Sec. V. 



II. MODEL HAMILTONIAN 

To begin we consider a homogeneous dilute gas of 
fermionic atoms with hyperfine spin / = 3/2 in a box 
with volume V . The second quantized Hamiltonian of 
the system is given by 



H = [ dr[5:^l(r)(-|^V^ 



-P\)i^\{'r) 



(1) 



+ ^ C/aiA2A3A4V'1i(i-)V'1,(i-)V'A3(i-)V'A4(i-)], 

Al,A2,A3,A4 

where Va (A = -3/2,-1/2,1/2,3/2) is the atomic field 
annihilation operator associated with atoms in the hy- 
perfine spin state |/ — 3/2, my — A). The summation in- 
dices in ([ij run through the values —3/2, —1/2, 1/2, 3/2. 
px is the bare atomic energy for spin state A. We will 
consider an effective quadratic Zeeman shift such that 
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P-i/2 = Pi/2 Ei-nd p_3/2 = P3/2 as the linear Zeeman shift 
can be gauged away. If we consider the s-wave scattering 
only, then the interaction between atoms can be charac- 
terized by the coefficients C/A1A2A3A4 which are obtained 
from the two-body interaction model U = g^Po + 32^2- 
Here, Pp is the projection operator which projects the 
pair into a total hyperfine spin F state, gi? is the inter- 
action strength in the total spin F channel, which, in 
the calculation, will be replaced in favor of the s-wave 
scattering length ap via the regularization procedure: 



1 



gp Awh^ap V ^ h^k"^ 



-E- 



For Fermi gases, there is no s-wave interaction in the odd 
total spin {F = 1, 3) channels, since these channels are 
forbidden by Pauli's exclusion principle. 



We expand the field operators with plane wave function 



'0A = X^k -^ke / and complete the spatial integral 



(2) to the Hamiltonian into momentum space: 



E'^AkAj, Ak + [c2(aL+p/3Lpi^k/^k' + Ml'+pi^Lp/3kak') + co(a[,^k/3k-p/3kak' + A' 

k,k',p 



k,A 



+ y ("L+pMk-pA^kak' + /3k'+p4-pi'k/3k' + aL+p^-p'^kak' + ML+p^Lp/^kA^k')] 

I 



(3) 



Here, cq = (50 + giM'^V, C2 = (52 - go)/2^ and 
^Ak ~ h^k'^ /2m + p\. For notational simplicity, we 
have used a, /3, /i and v to denote the annihilation op- 
erators for the component with spin quantum number 
TO/ = 1/2, —1/2, 3/2 and —3/2, respectively. From 
the Hamiltonian above, we can see clearly that the terms 
with coefficients C2 describe spin-mixing processes such 
as that depicted in Fig. [TJ Thus spin mixing requires 
C2 7^ or (70 7^ ff2. In the case go = 52, spin mixing 
does not exist and the population in each spin compo- 
nent is individually conserved. Under such a condition, 
the interaction obeys SU{'i) symmetry [Til. IT^. 

As our focus is on the effect of spin-mixing interaction, 
we will consider the case with go ^ gi- Great simplifi- 
cation can be further achieved by assuming 32 = and 
(?o 7^ 0, in which case the second line of Eq. ([3]) vanishes. 
This is the case we will consider in this work. We note 
that the essential physics does not change qualitatively 
if52y^0. 



written as 

B = E(^^k + /iA)AtAk 

k,A 

-t- ^[A*;3_kak - AV_k/Xk + h.c] , (4) 

k 

where hx = go Ek(A[Ak)/2y and A = - C3). To 

diagonalize this Hamiltonian, we can perform the Bogoli- 
ubov transformation 

"k = Ukflk + Vkb\, /3_k = Mk&-k - fkflk 

Mk = SkUk + ikwlk' '^-k = SkW-k - ik^k (5) 

Here, the new coefficients should meet the following con- 
ditions to ensure the anticommutativity of new operators. 

|skP + |<k|'-l, and |7/kP + kk|'-l. (6) 

Then we can write down the Bogoliubov-de Gennes 
(BdG) equations as 

f Sale + ha- p. A Wuk\_p/uk\ 

V A* -£:^k -hi,+ ~^ J {v^ J -"^'{vi, J 



III. GROUND STATE PROPERTIES 



Throughout this work, we take temperature to be zero. 
We will first consider the mean-field ground state prop- 
erty of the system. In the case of g2 = 0, it is obvious 
from Hamiltonian ([3]) that superfluid pairing can occur 
between spin components (§,— f) and (^,— ^). If we 
denote Ci = Z]k(/^-kak) and C3 = X^ki'^-kMk), the 
Hamiltonian under the mean-field approximation can be 



/ ^Mk + - A -A \ f Sk\ f Si, \ , , 

1^ -A* -f.k-/i. + M A^kj " ' ^v^k 

where we have introduced the chemical potential fl into 
the equations. We assume that there are 2N particles 
in total, and the population in opposite spin states are 
equal, i.e., Na — Np and — N^,. Then we have N — 
X^kd'^^kP + l^kP) and the order parameter is given by 

^ = -^E("k<-^kt;;). (8) 
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FIG. 2: (Color Online) (a), (b) The momentum distribu- 
tion in spin components ±1/2 and ±3/2 for quadratic Zeeman 
shift p = (solid blue), p = 1 (dot-dashed orange) and p = 10 
(dashed red), in units of the Fermi energy Ef- In the calcula- 
tion we set ao = —4 for (al) and (a2), and ao = 1 for (bl) and 
(b2), where the dimensionless interaction strength is defined 
as ao = IQaokp/TT with kp being the Fermi wavenumber. (c) 
Order parameter A as a function of p. 



Without loss of generality, we take the effective 
quadratic Zeeman shift as Pa — = P±i/2 — and 
Pp. = Pv = P±3/2 = P ^ 0- For a given p and the s-wave 
scattering length oq, we solve the BdG equations self- 
consistently, and typical results are presented in Fig. 

The upper panels of Fig. [5] display the momentum dis- 
tribution in different spin components at various values of 
p. As p increases, the population in the spin components 
nif = ±3/2 decreases, and that in nif — ±1/2 increases. 
This can be easily understood from the energetic point 
of view. However, somewhat surprisingly, for large p, the 
momentum distribution in m/ — ±1/2 components ap- 
proaches a step function, exemplifying a normal Fermi 
sea. Indeed, Fig. ^c) shows that the magnitude of the 
order parameter decreases as p increases. Intuitively, one 
might have thought that for very large p, the population 
in the irif — ±3/2 components becomes negligible and 
the system reduces into a two-component spin-1/2 Fermi 
gas, which should exhibits a superfluid ground state at 
zero temperature. Results presented in Fig. [2] apparently 
contradicts this intuitive picture. 
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FIG. 3: (Color Online) The bound state energy of the Cooper 
pair as a function of p for dimensionless interaction strength 
(defined in Fig. [2]) ao = —4 (solid blue) and ao = 1 (dashed 
red), respectively. 



are given by 



i?^„^(k) = i?„^VoMk) + (v>„^(k') - i^^A^)) 



2V 

90 



where E^p = ^^-2Ep and E^,, ^ >±^-2Ef + 2p, and 
■00^ and ifi^v are pairing amplitudes between m/ = ±1/2, 
•nif = ±3/2, respectively. The presence of the Fermi 
sea imposes the restriction that k and k' lie outside of 
the Fermi sea, i.e., fc,fc' > kp- A negative value of the 
eigenenergy E means that the two particles form a bound 
pairing with binding energy \E\. The larger \E\ is, the 
more strongly the pair binds. 

After some mathematical manipulation, the two eigen- 
value equations (jH]) leads to the following single equation 
for E: 



1 



go 

2V 



E 



1 



1 



E — Eap E — El 



(10) 



which can be solved numerically. The solution of i? as a 
function oip is displayed in Fig. [31 from which one can see 
that, as p increases, \E\ tends to zero. In other words, 
the pairs becomes less and less bound. This result is 
consistent with the many-body result illustrated in Fig. [5] 
We therefore reach the conclusion that the energy mis- 
match induced by the quadratic Zeeman shift p, together 
with the spin-exchange interaction, tends to break the 
pair apart. This phenomenon is reminiscent of the ef- 
fect by a magnetic impurity on a spin-1/2 superconduc- 
tor [20l.[2l|. In this latter system, the magnetic impurity 
induces an energy difference between the two pairing par- 
ticles and has the tendency of destroying the pairing. 



IV. SPIN-MIXING DYNAMICS 



To gain some insights into this problem, let us consider 
the so-called Cooper problem where two extra particles 
with attractive interaction lie on the surface of a filled 
Fermi sea noninteracting atoms. The eigenvalue equa- 
tions describing the two interacting particles in k-space 



So far we have focused on the ground state of the sys- 
tem. Now, let us turn to the spin-mixing dynamics. Be- 
fore dealing with the many-body situation, it may be 
helpful to investigate the Cooper problem first. We take 
the initial state to be the ground state of the Cooper 
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system under an effective quadratic Zeeman field with 
p = \OEp. At t = 0, we suddenly turn this field off so 
that p — Q and the system starts to evolve. The dynam- 
ics of system is governed by Eq. © after replacing E on 
the left hand side with ihd/dt. The equations can be 
simplified if we redefine two quantities as follows: 

■0± (k, t ) = (k, t) ± (k, t)] , 

and the dynamical equations for p = can be rewritten 
as 



at m 



(11) 



m 



zfi^V-(k,0 = ^^-(k,t) + ^^^_(k',t)(12 



V 



Hence the two equations for ■(/'± are decoupled and the 
interaction term only appears in the equation for ip- . 

The above equations can be easily solved and results 
are presented in Fig. 21 Initially, due to the presence 
of the large quadratic Zeeman shift, almost all the pop- 
ulations are in spin states nif = ±1/2. Spin mixing 
dynamics is initiated by quenching this Zeeman shift at 
i = 0. We compare the cases for two different values of 
the interaction strength: a weak interaction with ao = —4 
and a strong interaction with oq = 1- In both scenarios, 
damping is observed in spin mixing, and stronger inter- 
action gives rise to a much faster damping. This can be 
intuitively understood as follows: The initial state can 
be regarded as a superposition of different eigenstates of 
the quenched Hamiltonian. The stronger the interaction, 
the larger the number of the eigenstates contained in the 
initial state. For i > 0, different eigenstates oscillate at 
different frequencies which results in the damping of the 
population dynamics. The more eigenstates are involved, 
the faster the damping. Therefore, such damping is a re- 
sult of the intrinsic multi-mode nature of the Fermi gas. 
For a system of spinor condensate near zero tempera- 
ture, as all the atoms occupy the same lowest-energy or- 
bitals, nearly undamped spin-mixing oscillations can be 
observed UsHH. 

We now turn our attention to the dynamics in the 
many-body setting. As in the Cooper problem, we pre- 
pare the system in the ground state with p = lOEp and 
at t = quench the quadratic Zeeman field to zero. The 
ensuing mean-field dynamics can be simulated by solving 
the time-dependent BdG equations, obtained by replac- 
ing the eigenenergies at the right hand side of Eqs. ([7]) 
with ihd/dt [111. At each time step, the order parameter 
will be updated as 
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k 



{u^{t)vUt) - su{t)tl{t)) . 



Representative results are shown in Fig. [5] The dy- 
namics exhibits qualitatively similar properties as in the 
Cooper problem: damping is observed in both the dyan- 
mics of the population and that of the order parameter. 
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FIG. 4: (Color Online) Populations in spin states mf = ±1/2 
and lUf = ±3/2 as functions of time (in units of h/Ep) at 
two different interaction strengths. 
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FIG. 5: (Color Online) Time dependence of population (up- 
per panels) and order parameter A (lower panels), in units 
of Ef, for two different interaction strengths. The horizon- 
tal lines in the lower panels represent the value of A in the 
ground state with p = Q. Time is in units of h/ Ef- 



and the stronger the interaction, the faster the damping. 
This can be understood using a similar intuitive argu- 
ment we presented for the Cooper problem. We note that 
such damping was observed in the experiment reported 
in Ref. [19]. Finally we remark that the ground state cor- 
responding to p = should have equal population in all 
four spin states and an order parameter value indicated 
by the horizontal line in the lower panels of Fig. [51 In 
the absence of any dissipation, as in our simulation, the 
system remains far away from the ground state. 



V. CONCLUSION 

To conclude, we have examined the spin mixing inter- 
action of a degenerate Fermi gas with four internal spin 
components. It is quite remarkable that in the presence of 
an effective quadratic Zeeman field that shifts relatively 
the bare energies of spin-(±l/2) and spin-(±3/2) states, 
the system becomes almost normal. This is analogous to 
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the effect of the Zeeman field that breaks the symmetry of 
the two spin components in a two-component (spin-1/2) 
Fermi gas. We have also investigated the spin-mixing dy- 
namics initiated by quenching the quadratic Zeeman field 
and show that, unlike in the case of a spinor condensate, 
damping will necessarily occm- in a many-body Fermi gas 
due to the intrinsic multi-mode nature of fermions. 
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